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Energy levels statistics following the Gaussian Symplectic Ensemble (GSE) of Random Matrix 
Theory have been predicted theoretically and observed numerically in numerous quantum chaotic 
systems. However in all these systems there has been one unifying feature: the combination of half- 
integer spin and time-reversal invariance. Here we provide an alternative mechanism for obtaining 
GSE statistics that is based on geometric symmetries of a quantum system which alleviates the need 
for spin. As an example, we construct a quantum graph with a particular discrete symmetry given 
by the quaternion group Q8. GSE statistics is then observed within one of its subspectra. 



In the 1950s and 1960s Wigner and Dyson pioneered 
the use of random matrices in modelling the statistical 
properties of the energy eigenvalues belonging to compli- 
cated quantum systems [J [5] . The techniques they de- 
veloped spawned a new field of mathematics which has 
since become known as Random Matrix Theory (RMT) 
and its application has spread far and wide to many ar- 
eas of Mathematics and Physics [3] . In particular it was 
later conjectured 0] that the high-lying quantum energy 
levels of classically chaotic systems are faithful to random 
matrix averages. 

One of the cornerstones of RMT is Dyson's three-fold 
way [2] , which groups quantum systems without geomet- 
ric symmetries into three distinct types. The first occurs 
if time-reversal invariance is broken, for example by a 
magnetic field, meaning the quantum Hamiltonian H is 
inherently complex. The remaining two appear if there is 
an antiunitary time-reversal operator T which leaves H 
invariant, i.e. [T,H] = 0. They are then distinguished 
by either = 1 or = — 1, in which case H is real 
symmetric or quaternion-real respectively. For chaotic 
systems, RMT makes predictions in all three instances by 
averaging over an ensemble of Hermitian matrices with 
the appropriate internal structure and Gaussian weighted 
elements. These are referred to as the Gaussian Unitary, 
Orthogonal and Symplectic ensembles (GUE, GOE and 
GSE) . We note that the number of symmetry classes can 
be extended to ten if additional anti-commuting symme- 
tries are present E] but this is beyond the scope of this 
letter. 

In systems without geometrical symmetries time- 
reversal invariance with = — 1, and hence GSE statis- 
tics, can only arise if the wavefunctions have an even 
number of components, commonly associated with half- 
integer spin. For such systems GSE statistics have been 
predicted and/or observed numerically in examples such 
as quantum billiards 0, maps [S] and quantum graphs 
[HI, and explained using periodic-orbit theory [TUl ITT] . 
However to date there has been no experimental obser- 
vation. 

For systems with geometric symmetries the situation 
becomes more involved. Here the Hilbert space decom- 
poses into subspaces invariant under symmetry transfor- 



mations, and the spectral statistics inside these subspaces 
depends both on the system's behaviour under time re- 
versal and on the nature of the subspace. For exam- 
ple 3-fold rotationally invariant chaotic quantum systems 
display GUE statistics within certain subspectra even if 
they are time-reversal invariant (T^HUj. 

What is remarkable, and has not been observed be- 
fore, is that Dyson's formalism also permits GSE statis- 
tics to be obtained without any need for spin, provided 
the quantum system has a certain kind of discrete sym- 
metry. More precisely, we observe GSE statistics within 
the subspectra of time-reversal invariant systems with 
7"2 = 1 that are associated to so-called pseudo-real ir- 
reducible representations of the corresponding symmetry 
group. Our specific system is a quantum graph with an 
internal symmetry given by the quaternion group Q8. 

Let us first recall some important concepts regarding 
discrete symmetries in quantum mechanics [15j . If g is 
a classical symmetry operation in, say, position space, 
a corresponding quantum mechanical operator U (g) can 
be defined by U{g)i{}{r) — i}{g~^{r)), and the Hamilto- 
nian of the symmetric quantum system commutes with 
this operator. The unitary operators U{g) form a rep- 
resentation of the symmetry group G, i.e. they satisfy 
U{g)U{g') = U{gg') for all g,g' € Q. It can then be 
shown that the operators U{g) have block-diagonal form 
in an appropriate basis of eigenf unctions of the Hamil- 
tonian. The blocks correspond to so-called irreducible 
representations (irreps) of the symmetry group. For dis- 
crete groups there are only a finite number of irreps which 
we label by a and the corresponding block size is the di- 
mension Sq of the irrep a. The Sa eigenf unctions corre- 
sponding to the same block are energy-degenerate. If we 
assemble them into an Sa-dimensional vector \a, n) then 
the operators U (g) act as 

U{g)\a,n} = M(^\gf\a,n} (1) 

where n labels different blocks belonging to the same ir- 
rep a and M*-"^ (g) is the matrix representing g in the 
irrep a. In this way the spectrum falls into subspectra 
associated to the different irreps. For example, in a sys- 
tem with mirror symmetry the symmetry group consists 
of the identity e and the reflection operator r. Then 
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there are two one-dimensional irreps with M(±)(e) = 1 
and M'^^)(r) = ±1 corresponding to wavefunctions even 
and odd under reflection. 

Like the behaviour under time reversal, all irreps may 
be classified into one of three types, depending on how 
they are related to their complex conjugate. Firstly if 
there does not exist a unitary matrix S such that 



(2) 



then a is said to be complex. Alternatively if ([2| holds 
and S = then a is real as all M'^°'\g) can simulta- 
neously be made real by some unitary transformation. 
Whereas if ^ holds and S = —S^ then an appropriate 
unitary transformation leads to a quaternion real form 



consisting of 2 x 2 blocks 



-6* 



with a, 6 e C, and 



the representation is called pseudo-real. 

This classification is important as it defines which spec- 
tral statistics appear within each subspace [2] . In partic- 
ular we will argue in the following that spectral statistics 
in subspectra associated to pseudo-real irreps in time- 
reversal invariant systems with — 1 are appropriately 
modelled by the GSE. The main point is to show that the 
time-reversal operator T is not the appropriate symmetry 
operator for the subspace, and instead one has to intro- 
duce a modified operator T that determines the symme- 
try properties of the subspace. For simplicity, we consider 
irreps with matrices consisting of a single quaternion-real 
2x2 block, and we take T as the complex conjugation 

operator. In this case (2 ) holds with S = J ^ Now, 

crucially, if \a, n) satisfies ([l]) then the state T\a, n) does 
not, since the M^°'\g) are not real for all g e Q. In 
contrast, using ([2|, we see that T\a,n) — ST\a,n) does 
satisfy ([l]). Moreover T also commutes with H because 
S serves to exchange the components of |a,7i), while H 
acts on each component individually. Hence T is the ap- 
propriate time-reversal operator to be considered in our 
subspace. 

However f satisfies f ^ = S'^T'^ = -T^ = -1, indicat- 
ing that our subspace belongs to the symmetry class as- 
sociated with the GSE. We also note that since = —1, 
the states ja, n) and T|q;, n) are linearly independent and 
have the same energy, leading to Kramer's degeneracy 

m- 

The present argument for the appearance of GSE 
statistics relies purely on identifying the appropriate 
RMT symmetry class, in the spirit of [4j. However 
one may extend semiclassical methods for systems with- 
out geometrical symmetries (see p!7HT5] and references 
therein) to explain why individual chaotic systems are 
faithful to these predictions [M] [20] . This builds upon 
earlier semiclassical work by Keating and Robbins [13] 
who incidentally predicted that subspectra associated to 
pseudo-real irreps show GOE behaviour. However they 
only investigated the so-called diagonal approximation. 



in which GOE is indistinguishable from GSE if one does 
not remove Kramer's degeneracy. 

To exemplify how pseudo-real subspaces can arise in 
systems without spin we turn to quantum graphs [21] 
[22] . Quantum graphs consist of vertices v connected by 
one-dimensional bonds b — {vi^V2). Each bond has a 
specified length L}, and the time-independent Schrodinger 
equation on each bond reads 



H^{x) 



dx"^ 



ijj{x) = Eijj{x), 



(3) 



where x defines the position on each bond. At the vertices 
the wavefunctions have to satisfy boundary conditions 
that make the Hamiltonian self-adjoint. For instance one 
can consider Neumann (or Kirchhoff) boundary condi- 
tions; these require that at each vertex the wavefunc- 
tions of all adjacent bonds are equal and their outward 
pointing derivatives sum to zero. 

We have chosen to use quantum graphs since their 
spectral statistics have been shown to agree with the 
corresponding random matrix predictions in the limit of 
large, sufficiently well-connected graphs [23l [24] (assum- 
ing that the bond lengths are rationally independent). 
Moreover, in practice numerical experiments agree well 
with RMT already for relatively small graphs [21) . 

All symmetry operations on a graph may be given in 
terms of permutations of the vertices. A permutation 
5 is a symmetry if it leaves the connectivity and the 
bond lengths of the graph invariant, i.e., if for every bond 
(ui, V2) there is a bond {gvi,gv2) and it has equal length. 

In order to observe GSE statistics we must choose a dis- 
crete group which admits a pseudo-real irrep, the smallest 
and thereby easiest to construct is the quaternion group 

Q8 := {±1, ±1, ± J, ±K ■.I^ = J^ = K^= UK = -1}. 

Here all group elements can be written as products in- 
volving two generators, for example / and J. One can 
show there are five irreps: Four real one-dimensional 
irreps given by M{I) = ±1, M{J) = ±1, and a 
fifth two-dimensional and pseudo-real irrep given by the 
quaternion-real matrices 



m(5)(/) 







and M(5)(J) 



1 

-1 



The aim is therefore to construct a quantum graph 
with this Q8 symmetry. Here we turn to a standard 
group-theoretical tool for visualising the structure of the 
group, known as the Cayley graph [2S]. Cayley graphs 
can be constructed for any discrete group by taking the 
group elements as vertices and connecting them by bonds 
related to the generators. (Multiple Cayley graphs can 
exist for the same group if there are different possible 
choices for the generators.) 

In the example of Q8, see Fig. [ija), we draw bonds 
between two vertices representing group elements if one 



3 



element can be obtained from the other by right multipli- 
cation with either / or J. When interpreting the result 
as a quantum graph we choose the same length Lj for all 
bonds related to / and similarly Lj for J. 

The resulting graph is symmetric w.r.t. left multipli- 
cation of all elements in Q8. For example, let us consider 
a bond of length L/ given by 6 = {g, gl) where g G Q8, 
then application of any element h G Q8 leads to the bond 
6' = {hg,hgl) = (g',g'I) which is also a bond in our 
graph with length Lj. 

However at present our graph is still too small to be 
well described by random matrix theory. A larger graph 
which still retains the overall symmetry can be obtained if 
we replace the vertices corresponding to group elements 
with subgraphs, see Fig. [l|b). These subgraphs must 
be identical for each group element, meaning they must 
have the same number of vertices, the same connections, 
and the same bond lengths for analogous connections. 
We also need bonds connecting different subgraphs (as- 
sociated to different group elements), and these have to 
be chosen symmetrically. To be specific let us label the 
vertices corresponding to the group element g by Vg^m, 
m G N. Then for each subgraph g at least one vertex 
Vg,m must be connected to a vertex Vgi^m' m the sub- 
graph gl. Due to symmetry the possible choices for m 
and m' and the corresponding bond lengths must be the 
same for all g. Further bonds, subject to the same sym- 
metry conditions, have to connect subgraphs g and gJ. 
We will discuss later that for larger subgraphs it is, in 
fact, advisable to connect the subgraphs by more than 
one bond in order to obtain a well-connected graph that 
displays RMT statistics. 





FIG. 1. (a) The Cayley graph of the quaternion group Q8 
with solid (black) bonds corresponding to the generator / and 
dashed (red) to J and (b) an example subgraph being placed 
at a vertex of the Cayley graph. 

Now, importantly, our graph has subspectra associated 
to each of the irreps of Q8. This includes the four ID 
irreps of degeneracy one and the remaining 2D, pseudo- 
real irrep of degeneracy four (two from the dimension of 
the representation and two from Kramer's degeneracy). 
The use of subgraphs serves to add complexity, meaning 
each subspectrum is expected to have RMT statistics and 
in particular the subspectrum associated to the pseudo- 



real irrep will have GSE statistics. 

The wavefunctions in our graph are scalar, allevi- 
ating any requirement for half-integer spin. However 
two-component wavefunctions (associated to the two- 
dimensional irrep) come into play if one attempts to iso- 
late the GSE subspectrum using a so called 'quotient 
graph' [251 US] • This is essentially a fundamental domain 
equipped with the correct boundary conditions, analo- 
gous to the splitting of even and odd functions in a sys- 
tem with a reflection symmetry by taking the half system 
with either Neumann or Dirichlet boundary conditions on 
the former symmetry line. 

We illustrate the construction of this quotient graph by 
starting from the Cayley graph with eight vertices, each 
representing one group element. The simplest way to 
form the quotient graph is to take an eighth of the graph 
containing one vertex and half of each of the four gener- 
ating bonds attached to it. We choose quite arbitrarily 
the vertex K and cut the bonds (—J, K), {K, J), (/, K) 
and {K,—I) in half at the points t;i, i;2, W3, W4 shown in 
Fig. [2] (a). These points are related by symmetry oper- 
ations, in particular the application of / takes the bond 
{K, J) to {—J, K) and hence the point V2 in the middle 
of the intervening bond to vi. Similarly the application 
of J takes (K, —I) to (/, K) and hence V4 to V3. 





FIG. 2. (a) Definition of the points vi,V2,V3,V4, and the fun- 
damental domain, (b) Identification of the points vi with «2 
and with V4 with vertex conditions as explained in the text 
leads to the quotient graph. 

I'll ^2, and U4 now form the boundaries of the quo- 
tient graph and we have to identify boundary conditions 
that isolate the subspectrum associated to the pseudo- 
real representation. To do so it is helpful to again com- 
bine the pairs of degenerate energy eigenfunctions as- 
sociated to the pseudo-real representation into vectors 
tpix) — {x\a,n), where x denotes a position anywhere 
on the graph. We now evaluate U{I)tj}{vi). The defi- 
nition of U{I) and the symmetries of the graph imply 
U{I)ip{vi) = ip{I^^vi) = ip{v2)- Combining this with 
U{I)ip{vi) = Af(-'^)(/)^^(t;i) (see Eq. jl])) we obtain 

^{V2) = M^''Hlf'>P{vi) ■ (4) 

A similar result holds for the first derivatives if we let 
the coordinates along the bonds increase in the directions 
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indicated by arrows in Fig. [2] In this case we obtain a 
relation as in Q also for points moved compared to vi 
and V2 by the same amount, and differentiating w.r.t. 
this amount yields 



V''(i;2) = M(5)(/)V(^;i). 



(5) 



Analogous reasoning for the points V3 and V4 gives the 
conditions 



tA'(«4) = A/(5)(J)V(«3) 



(6) 
(7) 



Hence we identify vi with V2 and with V4 up to multi- 
plication of the ^ with a matrix, see Fig. ^h). The rela- 
tions Q to ([7]) have the effect of isolating the pseudo-real 
representation, now with two-component eigenfunctions 
supported on an eighth of the original graph. The same 
ideas can also be applied if, as described above, the ver- 
tices of the Cayley graph are replaced by subgraphs to 
increase complexity and hence generate random matrix 
statistics. In this case the quotient graph displays GSE 
statistics. Similar 'quotient systems' can be constructed 
to isolate the subspectra of any system with discrete sym- 
metries. 

The conditions Q to ([t]) provide an additional, intu- 
itive argument for the agreement with the GSE: A wave 
propagating through the graph will pick up factors as- 
sociated to / and J at the vertices vi, V2, W3, and v^. 
Repeated traversals of these vertices allow to sample the 
whole Q8 group, a discrete subgroup of SU (2) which de- 
scribes spin. This is similar to a discrete spin mechanism 
which, as was first noted in the context of the Dirac op- 
erator in [TU], is capable of generating GSE statistics. 





FIG. 3. (a) Quotient graph containing two bonds with M(I) 
conditions and two with M (J) and (b) Nearest neighbour 
spacing distribution, averaged over 10 graphs with 10,000 en- 
ergy levels each. 

For a numerical check we calculated the spectrum as- 
sociated to the pseudo-real representation in the quotient 
graph displayed in Fig. [s] (a). We took an average over 
ten random realisations with bond lengths distributed 
uniformly between and 1. Fig. |3] (b) shows a good 
agreement with Wigner's GSE prediction [T] [inj for the 
distribution P{s) of spacings s between neighbouring en- 
ergy levels (normalized to yield an average spacing of 1). 



The distribution of each individual realisation differs only 
slightly from the mean. The choice of two bonds with / 
and J conditions corresponds to a better connectivity in 
the full QS-symmetric graph (as mentioned earlier) than 
with only one bond. We also investigated larger graphs 
and found that we obtain better agreement with RMT if 
they are sufficiently well-connected. This is similar as in 
the case of non-symmetric graphs. 

In summary, we gave theoretical arguments and pro- 
vided numerical evidence that GSE statistics can be 
observed in systems without spin if they have certain 
symmetry groups that allow for pseudo-real representa- 
tions, as for example the quaternion group Q8. Quantum 
graphs as proposed here and e.g. built using optical fibres 
or coaxial cables [27] could lead to a first experimental 
observation of GSE statistics, avoiding the requirement of 
half-integer spin. It would be interesting to identify other 
experimental realisations of symmetries with pseudo-real 
representations and also investigate geometrical symme- 
tries within the framework of the new symmetry classes 

The authors would like to thank R. Band, J. Harrison 
and U. Smilansky for helpful advice. 



[4 

[5] 
[6] 

[7 

[8 

[9 

[10 

[11 
[12 

[13 

[14 

[15 

[16 

[17 
[18 

[19] 



E. Wigner, Ann. of Math. 62, 548564 (1955) 

F. J. Dyson, J. Math. Phys. 3, 1199 (1962) 

The Oxford Handbook of Random Matrix Theory, edited 
by G. Akemann, J. Baik, and P. D. Francesco (Oxford 
University Press, 2011) 

O. Bohigas, M. J. Giannoni, and C. Schmit, Phys. Rev. 
Lett. 52, 1 (1984) 

M. Zirnbauer, J. Math. Phys. 37, 4986 (1996) 

A. Ahland and M. Zirnbauer, Phys. Rev. B 55, 1142 

(1997) 

M. V. Berry and R. J. Mondragon, Proc. R. Soc. Loud. 
A 412, 53 (1987) 

S. Keppeler, J. Marklof, and F. Mezzadri, Nonlinearity 
14, 719 (2001) 

J. Bolte and J. Harrison, J. Phys. A: Math. Gen. 36, 2747 
(2003) 

J. Bolte and J. Harrison, J. Phys. A: Math. Gen. 36, 
L433 (2003) 

P. Braun, J. Phys. A: Math. Theor. 45, 045102 (2012) 
F. Leyvraz, C. Schmit, and T. H. Seligman, J. Phys. A: 
Math. Gen. 29, L575 (1996) 

J. P. Keating and J. M. Robbins, J. Phys. A: Math. Gen. 
30, L177 (1997) 

C. H. Joyner, S. Miiller, and M. Sieber, J. Phys. A: Math. 
Theor. 45, 205102 (2012) 

J. P. Elliot and P. G. Dawber, Symmetry m Physics, Vol- 
ume 1 (Macmillan, Basingstoke, 1979) 
F. Haake, Quantum Signatures of Chaos, 3rd ed. 
(Springer, Heidelberg, 2010) 

M. V. Berry Proc. R. Soc. Loud. A 400, 229 (1985) 
M. Sieber and K. Richter, Physica Scripta T90, 128 
(2001) 

S. Miiller, S. Heusler, A. Altland, P. Braun, and 



5 



F. Haake, New J. Phys. 11, 103025 (2009) 
[20] C. H. Joyner, Ph.D. thesis, University of Bristol (2012) 
[21] T. Kottos and U. Smilansky, Ann. Phys., NY 274, 76 

(1999) 

[22] S. Gnutzmann and U. Smilansky, Advances in Physics 
55, 527 (2006) 

[23] S. Gnutzmann and A. Altland, Phys. Rev. Lett. 93, 
194101 (2004) 

[24] G. Berkolaiko, in "Quantum Graphs and Their Appli- 
cations", Contemp. Math. No. 415, edited by S. F. 



G. Berkolaiko, R. Carlson and P. Kuchment (AMS, Prov- 
idence, RI, 2006) 

[25] R. Band, O. Parzanchevski, and G. Ben-Shach, J. Phys. 
A: Math. Theor. 42, 175202 (2009) 

[26] O. Parzanchevski and R. Band, J. Geom. Anal. 20, 
439471 (2010) 

[27] O. Hul, S. Bauch, P. Pakonski, Savytskyy, 
K. Zyczkowski, and L. Sirko, |Phys. Rev. E| 69, 056205 
(May 2004) 



